orci <- function(a, b, c, d, alpha) {
  # Calculate Odds Ratio 95% CI
  if(nrow(dt)==2 && ncol(dt)==2){
    a <- dt[1, 1]
    b <- dt[1, 2]
    c <- dt[2, 1]
    d <- dt[2, 2]
  
    or_est <- (a*d) / (b*c)
    tmp <- qnorm(1-alpha/2) * sqrt(1/a + 1/b + 1/c + 1/d)
    c1 <- exp(log(or_est) - tmp)
    c2 <- exp(log(or_est) + tmp)
  
    return(c(c1, c2))
  }else{
    return(0)
  }
    
}
